Self-Propulsion of Droplets by Spatially- Varying Roughness 
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Under partial wetting conditions, making a substrate uniformly rougher enhances the wetting 
characteristics of the corresponding smooth substrate - hydrophilic systems become even more hy- 
drophilic and hydrophobic systems even more hydrophobic. Here we show that spatial texturing of 
the roughness may lead to spontaneous propulsion of droplets. Individual droplets are driven toward 
regions of maximal roughness for intrinsically hydrophilic systems and toward regions of minimal 
roughness for intrinsically hydrophobic systems. Spatial texturing can be achieved by wrinkling the 
substrate with sinusoidal grooves whose wavelength varies in one direction (inhomogeneous wrin- 
kling) or lithographically etching a radial pattern of fractal (Koch curve) grooves on the substrate. 
Richer energy landscapes for droplet trajectories can be designed by combining roughness texturing 
with chemical or material patterning of the substrate. 

PACS numbers: 47.55.dr, 68.08. Be 
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Consider a liquid droplet partially wetting a solid 
substrate such as glass in contact with a gas such as 
air. Broadly speaking a substrate may wet easily (hy- 
drophilic) or poorly (hydrophobic) depending on the na- 
ture of the substrate, the liquid and the gas. More specifi- 
cally the three relevant interfacial surface tensions deter- 
mine the contact angle made by the liquid-gas contact 
line meeting the plane of the substrate. The contact an- 
gle is less than 90° for hydrophilic systems and greater 
than 90° for hydrophobic systems. A totally wetting thin 
film corresponds to vanishing contact angle and a com- 
plete spherical drop balanced at a point on the substrate 
corresponds to the super hydrophobic limit with a 180° 
contact angle. 

Uniform surface roughness amplifies the basic wetting 
characteristics of the corresponding planar system. For 
hydrophilic/hydrophobic systems the greater substrate 
area, for a given planar projection, available on the 
rough substrate makes the wetting more/less favorable 
and lowers/raises the contact angle. What about sur- 
faces with spatially inhomogeneous properties? Although 
variable chemical patterning [1] and Lcidcnfrost droplets 
contacting hot surfaces with asymmetric sawtooth pat- 
terns @, El have been thoroughly explored, this letter 
addresses the energetic driving forces acting on droplets 
on a substrate with pure inhomogeneous roughness and 
no other variability. Droplets will spontaneously move 
around in the landscape of the surface topography max- 
imizing or minimizing the roughness for naturally hy- 
drophilic/hydrophobic systems respectively. Thus self- 
propelled droplets can be engineered to follow prescribed 
paths without external drive by appropriately designing 
the surface topography. 

The principle of driving liquid droplets via roughness 
gradients is simple. The free energy of a liquid droplet 
on a substrate is F = —IA C + agyA t + cr^yA^v, where 
I = o'sv — o~sl is the imbibition parameter [4j, &sl, 
o~sv and o~lv ar e the respective surface tensions between 
the three phases (Solid/Liquid/Vapor), A c is the contact 
area between the droplet and substrate, At is the (con- 
stant) total area of a substrate and Ar,y is the area of the 




FIG. 1: (a) A liquid droplet sitting on a rough substrate has 
more contact area with the substrate (A c > Ao) than the 
same droplet on an otherwise identical flat substrate, (b) A 
liquid droplet partially wetting a hydrophilic substrate with 
inhomogeneous roughness has a smaller contact angle at the 
rougher end (9 a) than at the smoother end (9b)- 



liquid-vapor interface, which is taken to be constant even 
when a droplet moves. The system of a liquid droplet on 
a rough substrate may also be viewed as a droplet on 
a fiat substrate with an effective imbibition parameter 
I e tf resulting from the roughness. I e tf is defined by 
IA C = I e ffAo, where Aq is the planar projection of the 
actual contact area. For rough surfaces, A c > Aq, as 
shown in FigQJa), and therefore I e ff/I > 1- Up to irrel- 
evant constants, the free energy of a droplet on a rough 
substrate is 



F = -Ief f (x)A Q , 



(1) 



where I e ff varies from place to place when the spatial 
roughness is inhomogeneous. The wetting characteris- 
tics of a substrate/liquid composite system determines 
the sign of / and therefore I e ff- A hydrophilic system 
is characterized by / > and an acute contact angle 
9 = arccos(//<7Ly). A hydrophobic system is character- 
ized by I < and an obtuse contact angle Q. When 
I, and so I e ff, is positive (negative) a substrate lowers 
(increases) its free energy when covered by a liquid. This 
spontaneously drives droplets on hydrophilic (hydropho- 
bic) substrates towards rougher (smoother) regions re- 
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FIG. 3: The ratio of the effective imbibition parameter to the 
physical imbibition parameter versus F = a\k 2 for a liquid 
droplet partially wetting a sinusoidally modulated substrate. 



FIG. 2: (a) A liquid droplet partially wetting a substrate with 
a uniaxial sinusoidally modulated roughness, (b) Schematic 
plot of sinusoidal grooves with wavenumber monotonically in- 
creasing in the direction orthogonal to the sinusoidal height 
profile. 



spectively. Eq.([T]) can also be used to understand the 
movement of droplets on a chemically heterogeneous sub- 
strates where I e f / (x) depends on the wetting character- 
istics of the chemical composition at the corresponding 
position on the substrate. 

The self-propulsion of liquid droplets on substrates 
with inhomogeneous roughness can also be understood in 
terms of the uneven distribution of the Laplace pressure 
across the droplet. The contact angle for a rough sub- 
strate (9 r ) is given by cos 9 r — r cos 9, with r = A c /A Q. 
Thus surface roughness amplifies the intrinsic wetting 
properties of the corresponding planar substrate. Take 
a droplet spanning a hydrophilic surface that is rougher 
on the left than on the right, as illustrated in Fig|T] (b). 
The contact angle is then smaller on the left than on 
the right: 9a < 9b- The mean curvature H at the B 
end thus exceeds that at the A end, leading to a Laplace 
over-pressure (P = 2aLvH) gradient from right to left 
driving the droplet to the rougher part of the surface. 
The reverse argument applies to a hydrophobic substrate, 
leading to motion towards the smoother part of the sub- 
strate. The self-propulsion of a droplet on a substrate 
with spatially varying roughness clearly requires the size 
of the contact disk between the droplet and the substrate 
to be larger than the typical size over which the rough- 
ness varies significantly. 

To be specific, consider a droplet on a uniaxial si- 
nusoidal substrate, as shown in FigUJa), realizable via 
wrinkled membranes @. The height of the substrate is 
represented by z{x) = afcCos(fca;), with translational in- 
variance along the y-direction. The roughness amplitude 
a*; will be taken much smaller than the maximum height 
of the droplet so that the shape of the liquid-air interface 
is unaffected by the shape of the substrate. For small 
amplitude roughness liquid droplets can be in complete 
contact with the substrate, since air pockets do not form 
underneath the liquid [4, y5J . The effective imbibition pa- 



rameter is given by 



I k J a v 



sin (ky), 



(2) 



where k = kR is the dimensionless wavenumber and 
r = a\.k 2 . Clearly T, arising from the gradient of the 
substrate height, is the parameter controlling the effec- 
tive imbibition. For R >> Eq.(2) simplifies to Q 



1 7T 



(3) 



where R(x) = f* /2 d9\fi~ + x sin 9. In this limit the 
effective imbibition parameter is dependent only on the 
product of the amplitude and the wavenumber of the 
sinusoidal substrate. Fig|3] is a numerical plot of the 
monotonic growth of I e ff/I vs. T. I e ff/I is doubled 
for r w 8. A simple gradient of the effective imbibi- 
tion parameter over the substrate can be achieved by 
varying the wavenumber k along the groove axis (y), as 
shown in the schematic FigfJ^b). Droplets will migrate 
to maximize/minimize the contact area for intrinsically 
hydrophilic/hydrophobic substrates. The magnitude of 
the driving force along the groove axis is proportional to 
the gradient of the effective imbibition parameter: 



V„J, 
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V 1 eff 



[R(T)-S(T)] 



din k(y) 
dy 



(4) 



where S(x) = 



tt/2 



and R(x) - S(x) 



IT X 

4 



— V O (x 3 ) . The driving force thus depends rather 

weakly on k{y) and vanishes as V approaches zero. 

Rough substrates may also be designed based on frac- 
tals [8( . Consider a substrate with etched grooves whose 
cross section is the lower half of the Koch curve, as il- 
lustrated in Figg] 9]. The Koch curve may be con- 
structed by starting with an equilateral triangle of side 
length ao and perimeter Lq = 3ao, then recursively 
adding equilateral triangles symmetrically on each line 
segment. After n steps, the perimeter of the new graph 
becomes L n = L (4/3) n with the length of each ele- 
mentary line segment being a n = a (l/3)'\ The Haus- 
dorff, or self-similarity, dimension of the Koch curve is 




FIG. 4: A substrate etched by fractal grooves. The cross- 
sectional shape of the grooves is the lower half of the Koch 
curve, as shown in the inset. 



du = In 4/ In 3 ~ 1.26 [T(|. Now consider a set of close- 
packed evenly aligned straight grooves constructed from 
an n-th order Koch curve. The distance between two 
neighboring grooves is twice the breadth of a groove. The 
contact area between a droplet of radius R and the fractal 
substrate, in the limit R much bigger than the breadth 
of a groove, is A c = ttR 2 [1 + q(n)], where q(n = 0) = 3/4 
and q(n > 0) = 2(4/3) n ~ 2 - 1/2. Since q(n) is always 
positive, A c /Aq > 1. The effective imbibition parameter 
is 



l eff 
I 



= l + q{n). 



(5) 



Thus the effective imbibition parameter depends only on 
the order of the Koch curve n and is independent of the 
seed side length ao . I e f / /I is two for n — 1 and exceeds 
10 for n = 8. Note that the extra volume of liquid inside 
the finer structure of the fractal grooves is negligible in 
the limit of large n as the area A n of the n-th order Koch 
curve converges: liuin^aoAn = 2-^00/^ ~ a o « a o^> 
where h is the droplet thickness. 

An effective roughness gradient can be made by etch- 
ing a radial array of grooves with Koch cross-section, as 
sketched in Figj4] We take an intrinsically hydropho- 
bic surface inclined at an angle a to the horizontal. 
A droplet sitting near the origin (O) of the radial ar- 
ray will move upward to reduce the contact energy pro- 
vided the gradient of the effective imbibition parame- 
ter is sufficient to overcome gravity. The number of 
grooves covered by a droplet of size R is N = 2R/ (L50), 
where 59 is the angular distance between neighboring 
grooves and L is the distance of the droplet from the 
origin. For simplicity, consider the "far-field" limit in 
which the droplet is sufficiently far from the origin that 
all the grooves under the droplet are effectively paral- 
lel. The change in contact area due to the roughness is 
AA = A c — Aq = ^NTr(b n — Iq)R, where b n is the area 
of a groove of unit length with b n (n > 0) = ^-(f )" and 
b n ( n = 0) = 5ao/3, and Iq = 2ao/3 is the groove breadth. 
This results in an effective imbibition parameter 



eff 



= 1 



b n — lo 
L69 



A droplet rising up a distance 5L increases its grav- 
itational potential energy by SW = mgSL sin a, where 
m ~ pR 2 h is the mass of the droplet. Meanwhile 
the surface energy decreases by 5F = A 5I e ff, where 
5I e ff = I e ff(L+5L) — I e ff(L). Spontaneous climb there- 
fore requires 5I e tf/5L > mgsiaa/Ao. Inserting Eq.([5|) 
for the effective imbibition parameter yields 



J° > 



mg sin a 



A Ilo 



(7) 



(6) 



This condition can be satisfied for large droplets near 
the origin on substrates with dense grooves. The right 
hand side of Eq.© is of order 10 sin a for L ~ cm, R ~ 
mm, J ~ mN/m, 59 ~ clq/L and h ~ O.li?. Since 
(b n — lo)/lo ~ 10 for n = 5, radially carved grooves made 
of the 5-th order Koch curves would generate sufficient 
roughness gradient to drive droplets uphill. 

There are several points in our analysis that may ulti- 
mately call for a more thorough treatment. Sharp sub- 
strate edges impede the motion of droplets via pinning 
of the triple line Adhesion hysteresis may also 

arise from the microscopic interactions between a droplet 
and the substrate [11]. These two effects are the main 
source of frictional energy dissipation [l3[ . We have ne- 
glected entirely viscous dissipation due to internal fluid 
flow within moving droplets [1, [l4| • 

Droplet flow driven by inhomogeneous surface rough- 
ness is strictly downhill according to the gradient of 
the height profile of the surface but one may also vary 
the chemical composition of the surface so that the in- 
trinsic surface tensions are spatially dependent. The 
combination of chemical and roughness patterning of- 
fers a rich variety of potential structures to obtain de- 
sired flow patterns. The African beetle Stenocara fog- 
basks by tilting forward into the early morning fog-laden 
wind of the Namib desert and collecting micron-sized 
water droplets on the smooth hydrophilic peaks of its 
fused overwings (elytra) Once a sufficiently mas- 

sive droplet is formed it rolls downhill against the wind 
to pool in textured hydrophobic waxy troughs and from 
there to the beetle's mouth. Surface structures modeled 
on the Stenocara wings have been synthesized by creating 
hydrophilic patterns on superhydrophobic surfaces with 
water/2-propanol solutions of a polyelectrolyte [1 61 ) . 
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